Nano— wires with surface disorder: 
Giant localization lengths and quantum— to— classical crossover 
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We investigate electronic quantum transport through nano-wires with one-sided surface rough- 
ness. A magnetic field perpendicular to the scattering region is shown to lead to exponentially 
diverging localization lengths in the quantum-to-classical crossover regime. This effect can be 
quantitatively accounted for by tunneling between the regular and the chaotic components of the 
underlying mixed classical phase space. 
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Transport through a disordered medium is a key issue 
in solid state physics which comprises countless applica- 
tions in (micro-) electronics and optics JJ. The ubiq- 
uitous presence of disorder plays a prominent role for 
the behavior of transport coefficients governing, e.g., the 
metal-insulator transition ^ . The interest in disordered 
media has recently witnessed a revival due to new exper- 
imental possibilities to study the 'mesoscopic' regime of 
transport where a quantum-to-classical crossover gives 
rise to a host of interesting phenomena . 

In most investigations a static disorder is assumed to 
be present in the bulk of a material. The strength and 
distribution of the disorder potential determine whether 
transport will be ballistic, diffusive, or suppressed in the 
localization regime In nano-devices the reduction 

of system sizes leads, however, to an increased surface- 
to-volume ratio, for which surface roughness can repre- 
sent the dominant source of disorder scattering. While 
random matrix theory (RMT) is successful in describing 
bulk disordered systems , its application to wires with 
surface disorder is not straightforward 0- 

In the present paper we study electronic quantum 
transport through a nano-wire in the presence of one- 
sided surface disorder and a magnetic field. We show 
both numerically and analytically that by increasing the 
number of open channels N in the wire, or equivalently, 
by increasing the wavenumber kp , the localization length 
^ increases exponentially. Using a numerical approach 
that allows to study extremely long wires we show an in- 
crease by a factor 10^ (Fig.0. Such a giant localization 
length falls outside the scope of RMT predictions, ^ cx iV, 
previously studied for this system 0|. Instead it can be 
understood in terms of the underlying mixed regular- 
chaotic classical motion in the wire. We find that the 
conductance through the wire is controlled by tunneling 
from the regular to the chaotic part of phase space. This 



process, often referred to as "dynamical tunneling" 0, 
has been actively studied in quantum chaos and plays an 
important role in the context of classically transporting 
phase space structures H, 0, 0, ^| ■ Here we estab- 
lish a direct quantitative link between the exponential 
increase of the localization length in mesoscopic systems 
and the suppression of tunneling from the regular to the 
chaotic part of phase space in the semiclassical limit. 

We consider a 2D wire with surface disorder to which 
two leads of width W are attached (Fig. ^ inset), with a 
homogeneous magnetic field B perpendicular to the wire 
present throughout the system. We simulate the disor- 
der by a random sequence of vertical steps. The wire can 
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FIG. 1: (Color online) Localization length ^ for a wire with 
surface roughness vs kpW/n = l/hctt- Results are compared 
for wires with (a) one-sided disorder (OSD) with B (red 
■ ), (b) OSD with B = (green •), and (c) two-sided disorder 
with B ^ (blue A). In (a) an exponential increase of ^ is 
observed in excellent agreement with Eq. which has no 
adjustable parameters (dashed line). 
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thus be assembled from rectangular elements, referred 
to in the following as modules, with equal width I, but 
random heights /i, uniformly distributed in the interval 
[W — 6/2, W + S/2]. This particular representation of 
disorder allows for an efhcient numerical computation of 
quantum transport for remarkably long wires L — > cxd 
by employing the modular recursive Green's function 
method We first calculate the Green's functions 

for M = 20 rectangular modules with different heights. A 
random sequence of these modules is connected by means 
of a matrix Dyson equation. Extremely long wires can 
be reached by implementing an "exponentiation" algo- 
rithm |l4j : Instead of connecting the modules individu- 
ally, we iteratively construct different generations of "su- 
permodules" , each consisting of a randomly permuted 
sequence of M modules of the previous generation. Re- 
peating this process leads to the construction of wires 
whose length increases exponentially with the number of 
generations [l^ . 

The transmission (i„m) and reflection amplitudes 
[fmn) for an electron injected from the left are evaluated 
by projecting the Green's function at the Fermi energy 
E-p onto all lead modes to, 71 e {1, ... , N} in the entrance 
and exit lead, respectively. Here N — [fcpW/TrJ is the 
number of open lead modes and kp the Fermi wave num- 
ber. We obtain the localization length ^ in a wire com- 
posed of L modules (i.e. length LI) by analyzing the di- 
mensionless conductance g = Tv{tH) in the regime g <^1, 
extracting ^ from (Ing) ^ —L/S^. The brackets (. . .) in- 
dicate the ensemble average over 20 different realizations 
of disorder and 3 neighboring values of wave numbers fcp. 

For increasing fcp we adjust the magnetic field B such 
that the cyclotron radius Tc — hkp/{eB) remains con- 
stant. This leaves the classical dynamics invariant and 
allows for probing the quantum-to-classical crossover as 
fcp — > 00. We choose Tc = 3W and a disorder amplitude 
S — {2/3)W such that we obtain a large regular region in 
phase space (see below) and use a module width I = W/5. 
We find for one-sided disorder an exponential increase of 
the localization length ^ (Fig.^l, while ^ remains almost 
constant when either (i) the magnetic field is switched 
off or (ii) a two-sided disorder is considered. The lat- 
ter clearly rules out that the observed giant localization 
length is due to edge states of the quantum Hall effect Q . 

Before giving an analytic determination of the expo- 
nentially increasing localization length, we provide an ex- 
planation invoking the mixed classical phase-space struc- 
ture which captures the essential features of this increase. 

The classical dynamics inside the disordered wire is 
displayed by a Poincare section in Fig. ^jp, for a vertical 
cut at the wire entrance {x — 0) with periodic boundary 
conditions in the x-direction. The resulting section (y, 
Py ) for px > shows a large regular region with invariant 
curves corresponding to skipping motion along the lower 
straight boundary of the wire. Close to the upper dis- 
ordered boundary (y > W — S/2) the motion appears to 
be chaotic for all Py. A corresponding Poincare section 
for px < (not shown) is globally chaotic. The low- 




FIG. 2: (Color online) (a) Nano-wire with the regular trans- 
verse modes (green) m = 4,3,2,1 for k^W/TT = 14.6. The 
gray shaded part indicates the j/-range affected by disorder, 
(b) Poincare section showing a large regular island with outer- 
most torus (dashed), a chaotic sea (blue dots) and quantized 
tori corresponding to the regular modes (green) . (c) Poincare- 
Husimi functions of these modes and their quantizing tori. 



est transverse modes (Fig.|2K) of the incoming scattering 
wave functions overlap primarily with the regular island 
(Fig. I^b) . Only their exponential tunneling tail through 
the diamagnetic potential barrier (in Landau gauge) 



(1) 



touches the upper disordered surface dX y > W — 5/2. 
In Eq. JQ), TOg is the electron mass, ojc the cyclotron 
frequency, and the guiding center coordinate. These 
regular modes can be semiclassically quantized as [isl IT^ 

A BA , , , 

— = -— = (to- 1/4) with TO = 1,2, ... , (2) 
h h/e 

where A is the area in the Poincare section enclosed by a 
quantized torus and A = r^A/p-p is the area in position 
space enclosed by a segment of a skipping orbit. One 

finds A — pfVc arccos(l — i^) — (1 — i^)-\/l — (1 — u)'^ 

for < < i^max < 1, where vvc is the y-position at 
the top of the cyclotron orbit. The size Aj-cg of the reg- 
ular island is found for i/ = t'max = {W — S/2)/rc. The 
Poincare-Husimi projections (i.e. projections onto coher- 
ent states of the transverse eigenfunctions show, indeed, 
a density concentration near the quantized tori residing 
in the regular region of phase space (Fig. |2t). 

The lowest mode to = 1 in the center of the island has 
the smallest tunneling rate d, 0, 0| 



71 ^ exp -C 



A 



rcg 



(3) 



to the chaotic sea with some constant C (see below). 
Its temporal decay exp(— 7ii) together with its veloc- 
ity vi = hkx/me lead to an exponential decay as a 
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function of propagation length x, exp(— jix/ui). This 
gives a localization length C ~ Il3- When in- 

creasing fcp, while keeping the cyclotron radius Tc fixed, 
the classical dynamics remains invariant while the island 



area scales as A^^g — Orog^ps- Here ^pg — 2ppW is 
the area of the Poincare section and Orcg is the rela- 
tive size of the island. This semiclassical limit is thus 
equivalent to decreasing the effective Planck's constant 
hcs := h/Aps — {k-pW/Tr)~^ and results in an exponen- 
tial increase of the localization length 



C - exp C 



(4) 



for hcff — > 0, qualitatively explaining Fig. More- 
over, this exponential increase should set in when the 
first mode fits into the island, i.e. for Areg/h « 1. For 
the parameters of Fig. ^ we have I'max = 2/9, result- 
ing in the critical value k-pW/i: « 3.5 which is in very 
good agreement with the numerical result. By contrast, 
for two-sided disorder or for _B = no regular island 
with skipping orbits exists and ^ shows no exponential 
increase, see Fig. ^ 

We now turn to an analytical derivation of the localiza- 
tion length using the specifics of the scattering geometry 
(Fig.[T]inset). To this end we first calculate the transmis- 
sion amplitude tn of the transverse regular mode m — 1 
by considering its consecutive projections from one mod- 
ule to the next 

L-i 

ill ~ n / XhU){y) Xhu+i)iy) '^y ^ (5) 

where Xh(j){y) is the mode wave function in module j 
with height h(j). Eq. Q amounts to a sequence of sud- 
den approximations for the transition amplitude between 
adjacent surface steps. As the wave function is expo- 
nentially suppressed at the upper boundary, the scale 
Ik-p introduced by the corners drops out of the calcula- 
tion. For simplicity, a few technical approximations have 
been invoked, whose accuracy can be checked numeri- 
cally: i) terms in the transmission from one module to the 
next that involve reflection coefhcients and are typically 
smaller by a factor of 5 are neglected, ii) contributions 
from direct coupling between different island modes are 
neglected, and iii) the factor (27/— y ^. — y ° / . ^ ^ ) from the 
orthonormality relation for the x functions |13| is omitted 
in the above integral as its contribution is negligible. 

The modes pertaining to different heights h can be 
written as Xhijj) = [Xoo{y) - eh{y)]/Nh, where Xoo(y) is 
the mode wave function if there was no upper boundary, 
Sfi (y) is the correction that is largest at the upper bound- 
ary (where Xhih)=0), and Nh is a normalization factor. 
Keeping only terms of order 0{eh) and using a WKB ap- 
proximation for £/i(y) around y — /imin = — (5/2 leads 
to 



tn « (1 - (7)2-^/*=^ with (J = 



(^min) 



(6) 



(9) 
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FIG. 3: (Color online) (a) Averaged conductance (g) vs length 
L of the wire. The step-wise decrease is accompanied by the 
disappearance of the regular modes and the flooding of the 
island region by chaotic states. The Poincare-Husimi distri- 
butions to the left (right) of the curve correspond to scatter- 
ing from left to right (backscattering from right to right), (b) 
Transmission (Tm) of the incoming mode m vs L. 



According to Eq. ^ the coupling strength is quantita- 
tively determined by the tunneling electron density at 
y = hmin in the classically forbidden region of the ID 
diamagnetic potential, Eq. (Q. The conductance in the 
regime g ^ 1 is now given by 

5= |iiiP«exp(-4aL/M), (7) 

resulting in a localization length ^ = M/ (4cr) . Using a 
WKB approximation for Xoo iy) we find 



a h 



-2/3 
cff 



b] exp 



cK,\l-dh%'rf') , (8) 



with coefficients a = {WTr^y/^CMrjp b — 

-27rzoCM, c = 7r(32/9)i/V/'p''/'[l + (3/20)?7p-i], 
d = -zop^/V(2i/'V2/37/). Here zq w -2.338 is the 
first zero of the Airy function Ai(z), ( = Ai(2:)^d2:, 
V = hu^/W, and p = Tc/W are dimensionless param- 
eters [23. Eq. is in very good quantitative agree- 
ment with the numerically determined localization length 
(Fig. ^) . We conclude that tunneling from the regular 
phase space island is primarily due to interaction of each 
regular mode with the rough surface rather than via suc- 
cessive transitions from inner to outer island modes. 

We note that the constant C in Eqs. and is 
found to be C = 27r[l -f- (289/960W"^], which differs 
from C = 27r ^'^d C = 3 — In 4 ^3 derived for other 
examples of dynamical tunneling from a resonance-free 
regular island to a chaotic sea. We also note that the 
scaling behavior of ^ in Eq. © is reminiscent of previ- 
ously obtained results for diffusive 2D systems (see Q). 

For the case of a constant magnetic field B, increasing 
k-p increases the cyclotron radius, Tc oc /cp, and the classi- 
cal dynamics is no longer invariant. In particular the area 
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of the regular island Arog ^ Vkp shrinks compared to 
Aps ~ fcp as skipping motion is increasingly suppressed. 
Nevertheless, the arguments leading to Eqs. Q and (jHl 
remain applicable and yield a localization length that in- 
creases dramatically as ^ ~ exp (const ^/kp) in agreement 
with numerical observations (not shown). 

Now we turn to the behavior of the conductance for 
wires of lengths smaller than the localization length. 
Modes with larger m have larger amplitudes near the 
rough surface and thus couple more strongly to the 
chaotic part of phase space. They have, consequently, 
larger tunneling rates 7^ and smaller localization lengths 
~ ■ The successive elimination of modes as a func- 
tion of the length L of the wire results in a sequence of 
plateaus (Fig. . For L > the mode m no longer 
contributes to transport, as can be seen by its individual 
contribution to the transmission in Fig. ^jp. This dis- 
appearance of regular modes is reflected in the averaged 
Poincare-Husimi distributions calculated from incoher- 
ent superpositions of all modes entering from the left and 
scattering to the right. Also shown are the complemen- 
tary distributions obtained for backscattering from right 
to right. For small L these Poincare-Husimi functions 
are outside the regular island, while with increasing L 
they begin to "flood" it This process is complete 



for lengths L ^ ^. The complementarity of the Husimi 
distributions illustrates that tunneling between the reg- 
ular island and the chaotic sea proceeds symmetrically 
in both directions, as required by the unitarity of the 
scattering matrix. 

Summarizing, we have presented a numerical compu- 
tation and an analytical derivation for the exponential 
increase of the localization length in a two-dimensional 
system of a quantum wire with one-sided surface disor- 
der. Our approach, based on a mixed phase-space anal- 
ysis, also explains the increase of ^ over one order of 
magnitude under increase of the magnetic field observed 
in Ref. |6|. It sets in for a magnetic field for which the 
regular island is large enough to accomodate at least one 
quantum mechanical mode. Clearly, the RMT result, 
^ cc which ignores the mixed phase-space structure, 
no longer applies. Instead, we find that the giant lo- 
calization length (Fig. ^ in this disordered mesoscopic 
device is determined by the tunneling from the regular 
to the chaotic region, the rate of which is exponentially 
suppressed in the semiclassical regime. 
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